#name : fun binary exact ci caluclation
# key : bin.exact.ci.calc
# contributor: Shuguang Sun
# --
${1:binary_exact_ci} <- function(N = 20, resp = 5, alpha = 0.05) {

  exactlo <- qbeta(alpha / 2, resp, N - resp + 1)
  exacthi <- qbeta(1 - alpha/2, resp + 1, N - resp)

  list(exactlo = exactlo,
       exacthi = exacthi)
}

$1($0)